Pair Density Wave correlations in the Kondo-Heisenberg Model 
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We show, using density matrix renormalization group calculations complemented by field theoretic 
arguments, that the spin gapped phase of the one dimensional Kondo-Heisenberg model exhibits 
quasi-long range superconducting correlations only at a non-zero momentum. The local correlations 
in this phase resemble those of the pair density wave state which was recently proposed to describe 
the phenomenology of the striped ordered high temperature superconductor La2-iBa :E Cu04, in 
which the spin, charge, and superconducting orders are strongly intertwined. 



Recent experiments in the high temperature supercon- 
ductor La2_ K Ba a ;Cu04 near doping x = 1/8 have re- 
vealed a dramatic layer decoupling effect in which anoma- 
lous mesoscopic 2D superconductivity persists well above 
the macroscopic 3D superconducting transition tempera- 
ture, T c . [IH3] Moreover, the superconductivity coexists 
with static stripe (charge and spin) order. It has been 
proposed that the anomalous superconducting properties 
are evidence of the existence of a novel type of supercon- 
ducting state, the pair-density wave (PDW). [3HE] 

The PDW is a state in which charge, spin and su- 
perconducting (SC) orders are intertwined in a spatially 
modulated fashion. The SC order has a wave vector 
Q which is the same as that of the spin density wave 
(SDW) and half of the ordering wave vector 2Q of the 
charge density wave (CDW). Its SC order is Larkin- 
Ovchinnikov-like, but without the magnetization of the 
latter. Although much is known about the properties 
of this state[5J E], there is, as yet, no fully satisfactory 
microscopic theory. 

In the context of Bardeen-Cooper-Schrieffer type 
mean-field theories, a PDW is only ever stable at strong 
coupling [7] (i.e. outside the regime in which such treat- 
ments are reliable). Slave-boson mean-field theories of 
the t — J model find that, although the PDW is quite 
competitive energetically, it (barely) loses to the uniform 
d-wave SC state[8j. While early numerical variational 
Monte Carlo studies of the t — J model found a regime 
in which the PDW appeared to be stable [9], more recent 
studies have found that it has slightly higher variational 
energy than the uniform d-wave state [TUllllj . 

In the present paper we study the superconducting 
correlations in the ID Kondo-Heisenberg model (KHM). 
This is the simplest model in which one can investigate 
the interplay between strong antiferromagnetic ordering 
tendencies, represented by a Heisenberg chain, and possi- 
ble superconducting and charge-density wave orders, de- 
rived from an itinerant electron band to which it is cou- 
pled. The ID character of the model permits us to em- 
ploy the powerful numerical density-matrix renormaliza- 
tion group (DMRG)p2] and analytic bosonization meth- 
ods to solve the problem, despite the strong interactions. 



On the down side, there are special features of ID physics, 
which may raise questions concerning the applicability of 
the results to higher dimensional situations. On the other 
hand, especially since the order we are investigating is 
unidirectional, and thus has an essentially ID geometry, 
it is plausible that the local structure of the correlations 
up to intermediate scales are dimension independent. 

The key finding from our DMRG studies is that, for 
the range of parameters considered here, the ID KHM ex- 
hibits a spin-gapped phase with quasi-long range (power- 
law) PDW correlations,?. e. superconducting correlations 
which oscillate with a period 2b where b is the lattice con- 
stant of the Heisenberg chain. At the same time the uni- 
form singlet superconducting correlations are small and 
apparently fall exponentially with distance. Since the 
same model exhibits substantial, although short-ranged 
correlated, antiferromagnetic tendencies with the same 
period, this state can clearly be identified as a fluctuat- 
ing version of the long-sought PDW. Note that the occur- 
rence of a spin-gap in the ID Kondo-Heisenberg model 
has been discussed insightfully in the literature[T31 Q3] 
and the possibility of an oscillatory superconducting or- 
der parameter was previously inferred on the basis of 
bosonization studies. 15-20 However, we believe that 
this is the first place in which the existence and character 
of this state has been derived from a microscopic model, 
and the nature of the correlations is elucidated. |Hj 

Model: The ID KHM is defined as a one dimensional 
electron gas (1DEG) coupled to a spin-| chain: 

H = iflDEG + #Hcis + Hk (1) 



where 



HiBEG = -t^2c f j(T C j+la +h.C. -fl^Tlj, (2) 
#Heis = Jh / y Sj " Sj+i, (3) 
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Here, cl- ff creates an electron with spin a at site j, Sj 
is the spin ~ operator of the spin chain, and s a = |r a 
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FIG. 1: (Color online.) The spin gap vs. the electron con- 
centration in the 1DEG. Jh = Jk = 2i. The error bars are a 
result of the extrapolation to the thermodynamic limit. (Rel- 
ative to the extrapolation error, the DMRG truncation error 
is negligible.) 
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FIG. 2: (Color online.) (a) The SC order parameters <f> and <j>b 
(see text) as a function of position in an L = 64 system with 
n = 0.875. (b) Measurements of 4>pdw, <t>ie, 4>c vs. position. 
The oscillatory behavior close to the right boundary is an edge 
effect. 







are Pauli matrices). 



In typical physical circumstances in which Kondo 
physics arises, one would expect Jh and Jk "C t. In 
this limit, the length scales characterizing the Kondo ef- 
fect are exponentially large, and hence not readily ac- 
cessible by any numerical method. We therefore use 
Jh ~ Jk ~ t. On the basis of the field-theoretic analysis 
(see below), we expect the character of the phases to sur- 
vives to small J/t. Moreover, the Jh ~ Jk ~ t regime 
is not necessarily unphysical; it can be derived from the 
U — > oo limit of a Hubbard model on the spin-chain, with 
chemical potential chosen so that there is one electron per 
site, and with hopping matrix element along the chain, 
tn = y/JnU/2, and hopping between the spin-chain and 
the 1DEG, t K = y/J K U/2. 

Numerical results: The model (nj) was solved using 
DMRG on finite lattices with L = 32 — 128 and open 
boundary conditions. Up to m — 1800 states were kept, 
giving DMRG truncation errors smaller than 10~ 6 . 

Fig. fTlshows the spin gap A s = Eq (1) — Eq (0), where 
Eo (S z ) is the ground state energy of a system with a z 
spin projection S z . The spin gap was extrapolated to the 
thermodynamic (L — > oo) limit. The results are shown 
for Jh = Jk = 2i, as a function of the concentration 
of electrons in the 1DEG, n. Due to the particle-hole 
symmetry of the model, it is sufficient to consider n < 1. 

Near n = 1 there is a sizable spin gap 1 14]: A s decreases 
away from n = 1. For n = 1 (not shown), the spin 
gap is A s f=a 0.8, but there is also a finite charge gap in 
the L — > oo limit. We henceforth focus on n = 0.875, 
for which A s is substantial. Since A s persists at lower 
densities, we expect the low-energy properties at smaller 
n to be similar, although the correlation length is larger. 

PDW correlations: The opening of a spin gap is ex- 
pected to lead to enhanced SC (as well as CDW) corre- 
lations. To study these correlations, we have applied a 
local pair field to the left boundary [22]: 



#pair = A(c 



it 24. 



C 14. C 2t) 



h.c. 



(5) 



where we fixed A = 0.5i. (We have checked explic- 
itly that the results do not depend on the size of A.) 



The superconducting response of the system was probed 
by measuring the following induced order parameters 
throughout the system [2"3"]: 
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where <fi(j) and 4>b{J) are, respectively, the expectation 
of the singlet pair creation operator on site j and on 
the bond from site j to site j + 1. Fig. [2k shows <f){j) 
and 4>b(J) m an L = 64 system. <fi(j) appears to decay 
very rapidly away from the left boundary. 4>b(J) decays 
much more slowly, and exhibits pronounced oscillations 
as a function of position with wavevector q — n/b, as it 
changes its sign between every consecutive bonds. Longer 
periods are also apparent in the figure. These oscillations 
clearly indicate that the dominant pairing correlations 
are at a non-zero momentum. 

Refs. [T5Hl"7[ fi"9[ [20] proposed, based on bosonization, 
that the spin-gapped phase of the KHM has dominant 
pairing correlations at a non-zero wavevector, described 
by a "composite" order parameter [2D] 



oo = (-i) j <E4_ 1ct k S ) c 



u j + lcr 



]Sj 



(7) 



In addition, 
uniform (q = 



PDW order should be accompanied by a 
0) "charge 4e" order parameter 24J: 



04e (j) 



C n c ]i c j+n c i+^i> 



(8) 



Fig [2)3 shows 0pdw (j) = (-1)' 7 <Pb (j), as well as (j) c (j) 
and (j)ie (j), as a function of position, on a logarithmic 
scale. The largest, and most slowly decaying, order pa- 
rameter is 0pdw (i); suggesting that the system is best 
described as a fluctuating PDW state. As expected, 
4>4e (j) and 4> c (j) are non-zero, but small. <f> c (j) is mod- 
ulated as a function of position, while </>4 e (j) is smooth. 

The wavevectors of the leading SC and CDW fluctua- 
tions can be determined by a Fourier analysis of the SC 
and CDW orders. [2S] Fig. [3] shows the absolute values of 
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FIG. 3: (Color online.) Absolute values of the Fourier trans- 
forms of (a) the charge [{n(j)}] and (b) the SC [4>b(J)] orders. 
L = 32-128. 
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the Fourier transforms of <Pb (j) and n (j) 
for system sizes between L = 32 and 128. The charge 
density exhibits a large peak at q = 2kp + x/b which 
grows as a function of system size, where 2k p = xn. 
There are also small sub-leading features at q = 2fcp and 
q = x/b. The main feature in the Fourier transform of 
4>b is a pronounced peak at q = x/b, with a sub-leading 
peak at q = 2kp. This shows unambiguously that the 
dominant order in this system is a PDW with q = x/b, 
accompanied by CDW correlations at q = 2kp + x/b. 

In order to elucidate further the nature of the micro- 
scopic correlations in the system, we perform another 
simulation in which both a pair field [Eq ^} and a Zee- 
man field, Hz = —hS : - =1 (h — 0,5t) are applied to the 
left boundary of the system. The induced charge, super- 
conducting, and magnetic ({Sj}) order parameters are 
shown near the middle of the L = 64 system in Fig. [4^,. 
The magnetic order oscillates at wavevector q — x/b (the 
same as the PDW) with an envelope that decays expo- 
nentially on longer length scales. 

Next, we would like to understand what determines 
the PDW wavevector. We performed another calcula- 
tion in which the spin chain is "diluted", i.e. there is 
one spin site for every two 1DEG sites. Eq. ^ is 
replaced by # Hc is = JhJ2j^2j • &2j+2, and similarly 

H K = JKE jt a S 2j[4ja( sa )<T<T'C2 j<7 '}. Fig. [J illustrates 
the results for (Sj}, (rij), and 4>b near Ore middle of 
an L = 128 system. (In order to maintain a large spin 
gap, n was taken to be 0.625.) Clearly, the PDW order 
changes sign across every spin site, indicating that the 
dominant PDW wavevector is again q — x/b, where now 
b = 2. Thus, the period of the PDW is tied to that of the 
local (fluctuating) magnetic ordering. The local correla- 
tions in Fig. Hk-b are a one dimensional version of the 
phenomenologically proposed "striped-superconducting" 
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FIG. 4: (Color online.) Order parameters in (a) an L = 64 
KHM chain with n = 0.875, and (b) an L = 128 "diluted" 
KHM chain (with one spin site for each two 1DEG sites) with 
n = 0.625. Circles: 1DEG hole density (1 - (n(j))), bond 
color/thickness: the bond-centered SC amplitude 4>b(J). The 
arrows show the spin density (S z (j)). 



state for La2_xBa :E Cu04 |3]. 

Continuum limit: Analytical progress can be made in 
the limit Jh-X ^ Jk, where we may first take the contin- 
uum limit of both the 1DEG and the spin chain. We use 
a description in terms of the Bosonic fields <p c ,ip s and 
(p 8 , representing charge/spin fluctuations in the 1DEG 
and spin chain, respectively [and the respective conju- 
gate fields 6 C , 8 S and 9 S \. The Hamiltonian densities of 
the 1DEG and the spin chain take the form [T3l IT91 120] 

■HiDEG = ^2 v a[-^(d x 9 a ) + ------ {d x lfi N 



2K„ 



^Hcis = ^ v. [(dJs) 2 + (d x <f s ) 2 } 



(9) 



where K c , K s , v c , v s and v s are, respectively, the charge 
and spin Luttinger parameters of the 1DEG, and the 
corresponding charge and spin velocities. The vari- 
ous bosonized fields satisfy the commutation relation 
[tp a (x) , d x 6 a (a/)] — iS (x — x 1 ), and similarly for S , ip s . 
We neglect marginally irrelevant contributions to -ffHcis- 
For a an incommensurate filling n of the 1DEG, 
only "forward scattering" terms in the spin channel 
can couple the 1DEG and the spin chain. Up to ir- 
relevant (backscattering) operators, the Kondo Hamil- 
tonian density is H K = ^[(dx^+f - {d x ip-) 2 ] + 
H in t US where ± = i(0 a ± 9 S ), cp ± = -J=(£ fl ± 



V2' 



tp B ), and "Hint = gg ^^ ) [cos{y/Ix(p_) + cos(V4xtp + )]. 
(a is a microscopic cutoff.) Under renormalization, 
cos(V47r6>_) cos(\/4xip + ) is marginally relevant [13] . while 
cos(v / 4~7r6'_) cos(^/4xip^) is irrelevant, since it contains 
the dual fields 9_ and tp_. The strong coupling phase 
has a spin gap, while the charge degree of freedom (p c 
remains decoupled and gapless. 

Correlations in the spin gapped phase: The form of 
the dominant (slowest decaying) correlations follows from 
the following considerations. A theorem by Yamanaka 
et. al. [35] guarantees the existence of a charge zero, 
momentum 2k* F — xn to t gapless excitation, where ntot is 



the total electron density in the system (counting both 



the 1DEG and the spin chain). Here, n to 



t, there- 



fore 2kp = 2kp + ir/b. Let us denote the operator that 
creates these excitations Oik* ■ Since there is a spin gap, 
Oik* is necessarily a spin singlet, i.e., a CDW operator. 

In addition, as long as there is no charge gap, the 
singlet "77— pairing operator" [TH [20] 0,, — ip-ftp-^ — 
2^ exp[i\/27r (8 C — </> c )]> also creates gapless excitations. 
(V'i.o- annihilate right/left moving electrons with spin 
a =t, i, respectively.) This operator has total momen- 
tum —2kp and charge 2e. Therefore, the "PDW op- 
erator" Opdw = O^Oik* also creates gapless excita- 
tions. Adding the quantum numbers carried by O v and 
Oik* j we see that Opdw carries charge 2e and momen- 
tum ir/b. This guarantees the existence of quasi-long 
range PDW correlations in the spin gapped phase. As 
usual, the correlations of Opdw (as well as those of Oik* ) 
fall off with a non-universal exponent, which depends 
on K c . The (zero momentum) Cooper pair operator 
is Osc = ip+t^P-l = J ^e i ^( e < :+ ' l '°\ Its correlations 
are short ranged, since (j> s = (4> + + 0_) /v2; in the spin 
gapped phase the field #_ is pinned, while its dual </>_ un- 
dergoes strong fluctuations, suppressing the correlations 
of Osc ■ Consequently, the leading superconducting cor- 
relations are for operators with non-zero momentum. 

Generically, any singlet operator that carries charge 
2e and momentum ir/b is expected to couple to Opdw, 
and therefore to have quasi-long range correlations. For 
example, both 0pdw and <p c defined above have the cor- 
rect quantum numbers, and therefore their correlations 
should fall off with the same exponent as that of Opdw- 
According to our numerical simulations, the spin gapped 
phase has strong PDW correlations, so it is best char- 
acterized by the </>pdw order parameter. The results in 
Fig. [3] are fully consistent with the field theoretic analy- 
sis above. In particular, the density profile shows a large 
peak at q = 2k* F which grows with system size, indicating 
slowly decaying fluctuations centered at that wavevector. 
The pairing correlations are strongly peaked at q = ir/b, 
with a subdominant peak (which does not grow with L) 
at q = 2kp, corresponding to the gapless 77 pairing mode. 

Discussion: The correlations in the spin gapped 
phase of the ID KHM are best described as a PDW 
phase, which is a (quasi-)condensate of non-zero cen- 
ter of mass momentum Cooper pairs. Locally, the cor- 
relations are strikingly similar to those of the PDW 
state recently proposed to describe the striped phase of 
La2_ a: Ba :E Cu04, which intertwines spin, charge, and den- 
sity orders. A study of a two-chain KHM found, instead, 
dominant uniform pairing correlations 27J. It remains an 
important question whether the PDW state survives in 
other multi-chain generalizations of the present model. 
Finally, the ID KHM can be viewed as a variation of 
the three-band copper-oxide model [5Sj, with strongly lo- 
calized spins on the Cu sites and a 1DEG representing 



doped holes on O sites. Therefore it seems plausible that 
such a model can exhibit a PDW phase as well. Whether 
it can be realized in the physically relevant parameter 
regime remains to be seen. 
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Supplementary information 

Throughout the paper, we have used various boundary perturbations (either pairing fields, Zeeman fields, or the 
boundary itself) to excite various sorts of order parameters. The decay of these order parameters into the bulk gives 
information on the leading response of the system. In order to verify that our results (in particular, the fact that the 
dominant pairing fluctuations are at a non-zero wavevector) are not induced by the external boundary perturbations, 
we have made the following additional checks: 

1. We have varied the strength of the bond-centered boundary pair field, -ff pa ir, m the range A = 0.1 — 0.5t; 

2. We tried to apply an on-site pair field on the first site, rather than a bond-centered pair field, in order to verify 
that the nature of the correlations is not sensitive to the form of the pair fields; and 

3. We set A = and calculated the equal-time pair-pair correlation function. 

In the this Supplementary Information, we describe these results of these calculations. In all cases, the doping of 
the ID electron gas was n = 1.125 particles per site, and the following parameters were used: Jn/t = Jn/t = 2. 



Dependence on the boundary pair field strength 



Fig. [5] shows the induced pair fields 



^ U) 



<cU 



in three calcalutions, in which a bond-centered pair field 






1 



C 3l C 3+lV> 



H, 



pair 



= a (4t c U 



- c )A + n + H - c - 



(10) 



(ii) 



was applied to the first bond with a varying strength: A/t = 0.1,0.2,0.5. The system sizes were L = 64 in all 
calculations. The results are qualitatively similar to each other. In particular, the fact that the dominant response is 
at a non-zero wavevector q — n/a does not depend on the value of A/t. Note that, between A/t = 0.1 and 0.2, the 
response is roughly linear in the applied field. Between A/t = 0.2 and 0.5, the response nearly saturates. In all cases, 
we found that the on-site order parameter <f> (j) (not shown) decays very rapidly. 




FIG. 5: Induced superconducting order parameters in L 
strength, as a function of position along the chain. 



64 systems with bond-centered boundary pair fields of varying 



On— site boundary field 

In Fig. pi we plot <f> (j) and <pB (j) for an L = 64 system in which the following site-centered pair field was applied: 



H pah = Ac{ t c\ i +H.c, 



(12) 



with A = 0.5i. The results are very similar to the case of a bond-centered boundary pair field: the on-site order 
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FIG. 6: Induced superconducting order parameters, <j> and 4>b (see Eq. |10[ , in an L = 64 system with a site-centered boundary 
pair fields of strength A/t = 0.5. 

parameter 4> (j) decays very rapidly, while a bond-centered order parameter 4>b (j) is induced. 4>b (j) oscillates at a 
wavevector q = n/a and its amplitude decays slowly. This shows that the induced order parameter far away from the 
edges is not sensitive to the details of the edge fields. 



Equal— time correlations 

Finally, we calculate the pair-pair correlation function 

C {x) = {4> B ([L/2 - x/2]) 4> B {[L/2 + x/2])) 



(13) 



in an L — 16 system. [•] represents rounding to the nearest integer from below. The results (Fig[7| clearly show that 
beyond the first few neighbors, C (x) oscillates with a period of a, i.e. the dominant superconducting correlations are 
pair density wave correlations, in agreement with measurements of the response to edge fields. 
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FIG. 7: The equal-time two-point superconducting correlation function (Eq. 13 1 as a function of distance x. 
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